Ballistic transmission through a graphene bilayer 
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We calculate the Fermi energy dependence of the (time-averaged) current and shot noise in an 
impurity- free carbon bilayer (length L <C width W), and compare with known results for a mono- 
layer. At the Dirac point of charge neutrality, the bilayer transmits as two independent monolayers 
in parallel: Both current and noise are resonant at twice the monolayer value, so that their ratio 
(the Fano factor) has the same 1/3 value as in a monolayer — and the same value as in a diffusive 
metal. The range of Fermi energies around the Dirac point within which this pseudo-diffusive result 
holds is smaller, however, in a bilayer than in a monolayer (by a factor l±/L, with l± the interlayer 
coupling length). 

PACS numbers: 73.50.Td, 73.23.-b, 73.23.Ad, 73.63.-b 



-a 
c 

o 
o 



> 

(N 
&\ 
O 

\o 

O 

-i— > 
c3 



i 

O 
o 



X 



I. INTRODUCTION 

Undoped graphene has no free electrons, so an infinite 
sample cannot conduct electricity. A finite sample can 
conduct, because electrons injected at one end can be 
transmitted a distance L to the other end via so-called 
evanescent modes. These are modes that decay oc e~ L / x 
with a penetration depth A bounded from above by the 
width W of the sample. For a wide and narrow sam- 
ple (W ^S> L), there are many evanescent modes that 
contribute appreciably to the conductance. Because the 
transmission of an electron via an evanescent mode is a 
stochastic event, the current fluctuates in time — even 
in the absence of any scattering by impurities or lattice 
defects. Tworzydlo et ali found that the shot noise pro- 
duced by the evanescent modes in an undoped carbon 
monolayer (of length L <C width W) is pseudo-diffusive: 
The Fano factor F — P/2eI (ratio of noise power P and 
time-averaged current I) has the same value F = 1/3 as 
in a diffusive metal (while F = 1 for independent current 
pulses) £ 

A carbon bilayer has an additional length scale, not 
present in the monolayer of Ref. [l|, namely the interlayer 
coupling length It is an order of magnitude larger 
than the interatomic distance d within the layeri^^ 



, Tiv 3tii 
l± = — = - — d 
t± 2i_L 



lid 



(1) 



(with v w 10 6 m/s, d w 1.4 A, and t\\ « 3eV, re- 
spectively the carrier velocity, interatomic distance, and 
nearest-neigbour hopping energy within a single layer, 
and t± w 0.4 eV the nearest-neighbour hopping energy 
between two layers^) . Since L is typically large compared 
to the two layers are strongly coupled. In this paper 
we investigate what is the effect of interlayer coupling on 
the average current and shot noise. 

The model and calculation are outlined in Sees. [TTl and 
Mil Our main conclusion, presented in Sec. IIV| is that 
an undoped graphene bilayer has the same current and 
noise as two monolayers in parallel. The Fano factor, 
therefore, still equals 1/3 when the Fermi level coincides 
with the Dirac point (at which conduction and valence 



bands touch). However, the interval AEp ~ hvl±/L 2 in 
Fermi energy around the Dirac point where this pseudo- 
diffusive result holds is much narrower, by a factor l±/L, 
in a bilayer than it is in a monolayer. 

Our results for the mean current /, and hence for the 
conductance in a ballistic system, agree with those of 
Csertij^ but differ from two other recent calculations in 
a (weakly) disordered system*^ (The shot noise was not 
considered in Refs. A ballistic system like ours 

was studied recently by Katsnelson^ with different re- 
sults for both conductance and shot noise. We discuss 
the origin of the difference in Sec. |V] We conclude by 
connecting with experiments— in Sec. IVII 



II. MODEL 

We use the same setup as in Refs.fllfiol. shown schemat- 
ically in Fig. [1] A sheet of ballistic graphene in the x — y 
plane contains a weakly doped strip of width W and 
length L, and heavily doped contact regions for x < 
and x > L. The doping is controlled by gate voltages, 
which induce a potential profile of the form 



U(x) 



-Uoo if x < or x > L, 
if < x < L. 



(2) 



We use an abrupt potential step for simplicity, justified 
by the fact that any smoothing of the step over a distance 
small compared to L becomes irrelevant near the Dirac 
point, when the Fermi wave length > L. 

While Refs. fi]03 considered a graphene monolayer, 
governed by the 2x2 Dirac Hamiltonian, here we take a 
bilayer with 4x4 Hamiltonian^^ 



ti 



H 



U v(p x +ip y ) 
v(p x - ipy) U 

t± U v(p x - ip y ) 

v(p x + ipy) U 

(3) 

with p = —ihd/dr the momentum operator. 
The Hamiltonian acts on a four-component spinor 
(^Ai , , ' ^^2) with amplitudes on the A and B 
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FIG. 1: Schematic of the graphene bilayer. Top panel: Two 

stacked honeycomb lattices of carbon atoms in a strip between FIG. 2: Energy spectrum @ of the graphene bilayer, accord- 
metal contacts. Bottom panel: Variation of the electrostatic ing to the Hamiltonian ([3]). 
potential across the strip. 
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sublattices of the first and second layer. Only nearest- 
neighbor hopping is taken into account, either from A to 
B sites within a layer or between different layers. (Sites 
from the same sublattice but on different layers are not 
directly adjacent.) The Hamiltonian ([3]) describes low- 
energy excitations near one of the two Dirac points in 
the Brillouin zone, where conduction and valence bands 
touch. The other Dirac point and the spin degree of free- 
dom contribute a four-fold degeneracy factor to current 
and noise power. 

We have taken the same electrostatic potential U in 
both layers. In general, the potentials will differ ) 13 i 14 
but to study the special physics of undoped graphene it 
is necessary that they are both tuned to the Dirac point 
of each layer. This can be achieved by separate top and 
bottom gates (not shown in Fig. [T]). 

For free electrons in bilayer graphene, the relation be- 
tween energy s and total momentum k = (k 2 + ky) 1 ^ 2 as 
described by this Hamiltonian consists of four hyperbo- 
las, defined by 



e = T±t ± ± ^\t\ + k 2 , 



(4a) 
(4b) 



plotted in Fig. [2j (For notational convenience, we use 
units such that hv = 1 in most equations.) 

We calculate the transmission matrix t through the 
graphene strip at the Fermi energy, and then obtain 
the conductance and noise power from the Landauer- 



G = Go Tr tt* , P = P Trtt t (l-ii t ), 
Tr«t(i _ tit) 



F = 



(5) 
(G) 



with G = 4e 2 //i, P = 2e\V\G , and V the voltage 
applied between the contact regions. The results depend 
on the degree of doping in the graphene strip (varied by 
varying Ep), but they become independent of the degree 
of doping of the contact regions if Uoo ^> t± . 



III. TRANSMISSION PROBABILITIES 

We calculate the transmission matrix by matching 
eigenstates of the Hamiltonian © at the two interfaces 
x = and x = L. This procedure is similar to a cal- 
culation of non-relativistic scattering by a rectangular 
barrier in a two-dimensional waveguide. There are two 
differences. Firstly, the Hamiltonian © is a first-order 
differential operator, and hence only the wavefunction 
and not its derivative is continuous at the interface. Sec- 
ondly, the spectrum contains both positive and negative 
energy eigenstates. 

The eigenstates of H for U = have been given in 
Ref. [H. They may be characterized as follows. For given 
energy e and transverse momentum k v , we define two 
longitudinal momenta 



kx± = \/{e ± htj 



H 2 



4"-L 



U2 

Ky. 



(7) 



The square root is taken with argument in the interval 
[0, 7r). Associated with each real k x+ there are two prop- 
agating modes, one left-going (f>^ , and one right-going 
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4>^+- Two more propagating modes 4>^ _ and _ are 
associated with each real k x ^. These eigenstates of H 
are given by 



4>?Ax,y) =N± 



4>i±(x,y) = N± 



V „i(k m ±x+k y y) 



(8a) 



V k x ± 4- iky 



( 



±k x ±±ik y I e i(-fc a>± x+fc y y) ) ^) 



with JV± = (4M / efc x ±) 2 a normalization constant such 
that each state carries unit current 



I = ev I dy 
Jo 



a x 
a x 



(9) 



in the positive or negative x-direction. 

For each k y we have two left-incident scattering states 
ip £ .± at energy s. In the region x < to the left of the 
strip they have the form 

1pe,± = 0?+U^,± + r +( £ ,K)4>e+U x ,+ + r -( s iky)<i>e+U OQ -, 



while to the right of the strip (x > L) one has 



(10) 
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FIG. 3: Total transmission probability Tr(tt^) as a function 
of e and k y for Uac = 50 tj_ and L = 50 Zj_ . Darkly shaded 
regions indicate high transmission. Grey dashed lines indi- 
cate the estimate (114[) for the occurrence of resonances in 
regions (a) and (b), while solid lines indicate the boundary 
between propagating and evanescent modes. Arrows point to 
the resonances of evanescent modes close to the Dirac point, 
responsible for the pseudo-diffusive transport. 



Tp E ,±=t±(e,k v )cl>K +Uaot+ +t±(e,ky)<t>? +Uao! _. (11) 



IV. RESULTS 



For e the form of the solution in the region x £ [0, L] 
is self-evidently a linear combination of the four solutions 
^e± ' ^e± ■ Care must however be taken in analytical work 
to use proper linear combinations of these modes that 
remain linearly independent exactly at e = (the Dirac 
point). (See Appendix pA~l for explicit formulas.) 

The four transmission amplitudes t± for given e and 
k y can be combined in the transmission matrix 

^^-{t-^ky) tZ(s,k y )J- ^ 

We consider a short and wide geometry L <C IT, in which 
the boundary conditions in the y-direction become irrel- 
evant. For simplicity, we take periodic boundary con- 
ditions, such that k y is quantized as ky in = 2im/W, 
n = 0,±1,±2, .... In the regime L -C W, \e\ < 
considered here, both the discreteness and the finiteness 
of the modes in the contact region can be ignored. As 
a consequence, the traces in Eqs. ([5]) and ([6]) may be 
replaced by integrals through the prescription 

W f°° 

Tr (tit)? - _ / dk y V [T„(E F , k y )) p , (13) 
where T± are the two eigenvalues of ttf . 



Fig. [3] contains a grey-scale plot of the total transmis- 
sion probability Tr (tit) as a function of k y and s. Darkly 
shaded regions indicate resonances of high transmission, 
similar to those found in Ref. Il5l . 

The location e res of resonances can be estimated by 
equating k x L/n to an integer n. This yields the curves 

41(M = ± \M + (^) 2 + ^ (i4) 

indicated in the figure by dashed lines. It is seen that 
good agreement is reached for \k y \ <C 1/L and again for 
\k y \ S> 1/L. For \k y L\ ~ 1 there is a cross over. In regions 

(c) and (d), demarkated by the curves Eresj the transmis- 
sion generally drops to zero, since in these regions the 
longitudinal momentum k x is imaginary. 

There is however a curious feature close to e,k y = 0. 
The resonance closest to the Dirac point behaves differ- 
ently from all the other resonances. When \k y \ is in- 
creased, it moves closer to the Dirac point rather than 
away from it, eventually crossing into regions (c) and (d) 
of evanescent modes. It is this resonance of evanescent 
modes that is responsible for the pseudo-diffusive trans- 
port at the Dirac point. 

At e = 0, the exact formula for the eigenvalues of ttf 



ky [t±/flV] 



FIG. 4: Solid curves: Transmission coefficients T± of the bi- 
layer according to Eq. (|15p at L — 50 ij_ ■ These coefficients 
are displaced copies of the monolayer result (dashed). 



in the Uoo — *• oo limit is 
T±(e = 0,k y ) 

k c = — In 
L 



cosh =f A; C )L 




(15) 
(16) 



In Fig. d] the two transmission coefficients T±(0,k y ) 
are compared to the single transmission coefficient 



T, 



monolayc 



T (0,k y ) = 1/cosh (k y L) of the monolayer 



1.10 



Details of the calculation may be found in Appendix [X] 
Since the two bilayer coefficients are displaced copies 
of the monolayer coefficient, any observable of the form 
A = Tr f(tt'), with / an arbitrary function is twice as 
large in a bilayer as it is in a monolayer. From Eqs. © 
and p3p we obtain 
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FIG. 5: Conductivity a (top) and Fano factor F (bottom) of 
the bilayer, as a function of the Fermi energy Ef measured 
from the Dirac point for Uoo = 50 1± and L = 50 1±. Abrupt 
features occur at Ef ~ e;S }(k y = 0) [vertical lines, given by 
Eq. (US)]. 



DEPENDENCE ON THE POTENTIAL IN 
THE CONTACT REGION 



^bilayer — 2G m0 nolaycr 
-^bilayer — 2_P m0 nolayer ~ 



2G W 
4e|V|G W 



3tt 



^bilayer *i 



monolayer 



1 

3' 



(17) 

(18) 
(19) 



Figure [5] contains plots of both the conductivity a = 
GL/W and Fano factor of the bilayer around the Dirac 
point. At energies associated with resonances at normal 
incidence, 



41(0) = ± 



ir 2 hv 



l j-n* + 0{l^/Lf 



(20) 



the conductivity and Fano factor show abrupt features. 
The width AE F = 2eiU = 2Tr 2 hvlj_/L 2 of the energy 
window between the resonances that straddle the Dirac 
point in the bilayer is smaller by a factor l±/L than in 
the monolayer. 



So far we have assumed that the potential Uoo in the 
contact region is large compared to the band splitting 
t± near the Dirac point of the graphene bilayer. We be- 
lieve that this is the appropriate regime to model a nor- 
mal metal contact to the graphene sheet, which couples 
equally well to the two sublattices on each layer. 

It is of interest to determine how large the ratio Uoo/t± 
should be to reach the contact-independent limit of the 
previous section. Note that for Uoo > t± there are two 
left-incident propagating modes in the leads for each e 
and k y . When Uoo becomes smaller than t± one of the 
two modes becomes evanescent, leading to an abrupt 
change in the conductivity and the Fano factor. This 
is evident in Fig. [SI For Uoo —t± > hv/L, the conductiv- 
ity and Fano factor have almost reached their Uoo — * oo 
limits. For Uoo ^ t± the conductivity is smaller and the 
Fano factor larger than when Uoo > tx- Both quanti- 
ties vanish when the Fermi momentum y/U ot±/'f in the 
contact region drops below h/L and the contact region 
is effectively depleted of carriers. 
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FIG. 6: Dependence of the conductivity and Fano factor at 
the Dirac point on the potential Uca in the contact region, for 
L — 100 lx- Thin horizontal lines indicate the values of Ref. 
Il2l . The values obtained in this paper correspond to a plateau 
reached for Uoo/tx ^ 1. 



These finite- [Too results can be used to make contact 
with the previous calculation of Katsnelson^ who found 
a conductivity a = Gq/2 and a Fano factor F = 1 — 2/tt 
at the Dirac point, in the regime hv/L <^ y/Uootx <C tx- 
These values are indicated in Fig. [S] by horizontal lines. 
The intersection point with our curves occurs at nearly 
the same value of Uoo/tx for both quantities. The inter- 
section point moves closer and closer to Uoo — as the 
sample length L is increased, but there is no clear plateau 
around the intersection point. Moreover, as shown in Ap- 
pendix B, the intersection point does not correspond to a 



minimum or maximum as a function of the Fermi energy 
so that these values would be difficult to extract from a 
measurement. 

We do believe that the results of Ref. [H describe the 
asymptotic limit L/lx — > 00 at Ep = 0. however, be- 
cause in this limit the width AEp ~ Tivlx/L 2 of the 
resonance at the Dirac point vanishes, it seems unob- 
servable. 

VI. CONCLUSION 

In conclusion, we have demonstrated that the pseudo- 
diffusive transport at the Dirac point, discovered in Ref. 
1 for a carbon monolayer, holds in a bilayer as well. All 
moments of the current fluctuations have the same rela- 
tion to the mean current as in a diffusive metal. In par- 
ticular, the Fano factor has the 1/3 value characteristic 
of diffusive transport, even though the bilayer is assumed 
to be free of impurities or lattice defects. 

Although we found that an undoped bilayer transmits 
as two undoped monolayers in parallel, the two systems 
behave very different away from charge neutrality. The 
resonance of evanescent modes around the Dirac point 
of zero Fermi energy has width AEp ~ hvlx/L 2 in a 
bilayer, which is smaller than the width in a monolayer 
by the ratio of the interlayer coupling length lx and the 
separation L of the metal contacts. 

Since lx ~ 1.5 nm, one would not be able to re- 
solve this resonance in the ^m-size samples of Ref. 
[ill . These experiments found no qualitative difference 
in the conductance-versus-gate-voltage dependence of 
monolayer and bilayer graphcnc, both showing a min- 
imum conductivity at the Dirac point of Go- Smaller 
junctions in the 10-100 nm range as are now being fabri- 
cated should make it possible to resolve the transmission 
resonance of evanescent modes predicted here, and to ob- 
serve the unusual pseudo-diffusive dynamics associated 
with it. 



Acknowledgments 

This research was supported by the Dutch Science 
Foundation NWO/FOM. We have benefitted from dis- 
cussions with E. McCann. 



I 

APPENDIX A: TRANSMISSION EIGENVALUES AT THE DIRAC POINT 

In this Appendix we give some detail of the calculation that leads to the transmission coefficients T±(e = 0, k y ) of 
Eq. (|15p. At the Dirac point and in the limit of large Uoo, the left-incident eigenstates of the Hamiltonian (J3J) are of 
the form 

! e ik v y + ( r ±£L + r ±£L) g-Jt/^z] X < 0, 

e^v[(4xi + 4x2)e k * x + (4x3+4xi)e- k **] < x < L, (Al) 
e ik yVe iU^{x-L) [- t ±£fl + t ±£Rj X > L, 
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with the definitions 





(A2) 



Xi 



Xi = 




,X3 = 




Xi = 








(A3) 



These eigenstates must be continuous at x 
with 



M = 



( 1 

-1 

-1 
1 







-1 

1 1 

-1 
1 

z 









1 



-iLt±z 
z 




and x — L, leading to an 8 x 8 system of linear equations Mb± = c± 

i 



l 

o 




1 




-iLt i 



°\ 






-z 
-I 
-1 
-z 



( r \\ 



b+ = 



r_ 
ci 

C4 



Tl 
1 
1 




Vo / 



(A4) 



We abbreviated z = e 
t 



By solving these equations, one finds the transmission matrix 



2% 



2 + (L/l ± y + 2cosh(2k y L) 
The eigenvalues of ttf are then given by Eq. (fTS!) 



(L/l± ~ 2i) cosh(k y L) 
— (L/l±) sinh(fc !/ L) 



(L/l±) smh(kyL) 
{L/l± + 2i) cosh(kyL) 



(A5) 



r 



APPENDIX B: FOUR-BAND VERSUS 
TWO-BAND HAMILTONIAN 

In this Appendix we verify that the difference in the 
results obtained here and in Ref. [H is not due to the 
different Hamiltonians used in these two calculations. 

In Refill^ the limit tj_ — * oo was taken at the beginning 
of the calculation, reducing the 4x4 Hamiltonian (J3j> to 
the effective 2x2 Hamiltonian 4 



Only the two lowest bands near the Dirac point are re- 
tained in iJ ff, as is appropriate for the regime Uoo <C t±. 



H P ff = — - — 



(p x - ip y ) 2 

t± V (Px + iPy) 2 



-U(x) 



1 
1 

(Bl) 



We have repeated the calculation of conductance and 
Fano factor using both Hamiltonians ([3]) and (IBlj) . for 
parameter values corresponding to the intersection point 
of Fig. [6l and find good agreement (see Fig. [7]) . 
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FIG. 7: Conductivity (top) and Fano factor (bottom) around 
the Dirac point, for L — 100 lx and Uoa = 0.2 tx- (These 
parameter values correspond to the intersection point of our 
curves with the prediction of Ref.0 in Fig. [6]) The solid lines 
were obtained using the four-band Hamiltonian ((3| , while the 
dashed lines were obtained from the two-band Hamiltonian 
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